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I. INTRODUCTION 



The "academic" wind tunnel at NPS is inadequate for most 
tasks, due to the excessive level of turbulence at high 
speed. This tunnel was designed for a potential speed of 200 
Knots in the test section, but this value has not been 
reached yet. This project was undertaken with the purpose of 
evaluating, by a different method from that used for 
designing, the performance of the tunnel in order to 
establish possible causes of the problem. 

The two counter - rotating fans could be the cause of 
the trouble since they were designed using a blade element 
theory, that does not take account of the three dimensional 
effects that are present at medium and high solidities. 
These effects cause changes in the apparent velocity vector, 
therefore the blades were analysed by a " New Blade Element 
Method for Calculating the Performance of High and 
Intermediate Solidity Axial Flow Fans " due to Borst 
[Ref. 1: p.l]. This method determines an induced angle of 
attack which changes the apparent velocity past the airfoil 
and it has been found to be a reliable measure of the three 
dimensional effects. The data required for this method comes 
from the two dimensional airfoil data, that is available in 
the current literature. This method was used in reference 1 
(p.10) assuming zero drag, and the predictions showed 
excellent agreement with measurements. The small difference 
between the results can be attributed to the assumption of 
zero drag. The performance of the first and second stage 
blades of the academic wind tunnel is evaluated assuming non 
zero drag, so it is expected that the results will- be 
realistic . 

A possible cause of the inability to reach the specified 
speed could be excessive losses around the tunnel circuit, 
principally at the first diffuser and the corners. The first 
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diffuser has a high angle of divergence which can provide a 
flow separation, further losses and turbulence. The 
evaluation of losses around the tunnel is made by means of 
the method given by RAE & POPE [Ref. 2:' p.87], which breaks 
the tunnel down into sections and calculates the losses for 
each one. The four corner losses are evaluated by three 
methods : 

1 considering one third of the losses due to skin 
friction and two third losses due to_ rotation; an 
empirical relation given by RAE & POPE [Ref. 2: p.90], 
was used to account for the different types of corner 
vane , 

2 considering empirical values found by testing and 

3 considering variation of resistance coefficient with 

t ap/ chord ratio of thin corner vanes given by PANKHURST 
HOLDER [Ref. 3: p.93]. 

The first diffuser, since it has variable divergence, 
was broken into four parts and the losses evaluated part by 
part through the method given by RAE & POPE [Ref. 2: p.89], 
this shows that the losses are strongly dependent on the 
diffuser angle. 

After all losses were evaluated , a relation between the 
air inflow velocity to the blades and the power output from 
the blades is found and the non-stall operation envelope is 
determined. The temperature and pressure assumed by LARSON 
in the original blade design [Ref. 4: p.5], are T= 100° F 

O 

p= 2,246 (lbf/ft" 1 ) and, therefore the values of kinematic 

O 

viscosity (j>) of 0.000166 (ft* 1 / sec) and sound velocity of 
1160 ft /sec, as presented in Appendix C, will be used for 
evaluation of the Reynolds and Mach numbers. 
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II. EVALUATION OF LOSSES 



A. INTRODUCTION 

As mentioned in chapter one, we will use the procedure 
outlined by RAE & POPE [Ref. 2: p.87], that consists of 
breaking the tunnel down into cylindrical sections, 
expanding sections, contracting sections and corners, in 
order to evaluate the tunnel losses. 

The sections are numbered in such way that those 
sections with the same characteristics, that is, those that 
use the same loss equation, are kept together. This 
facilitates the evaluation by computer. Then, the 
cylindrical sections are numbered from 1 to 4, the expanding 
sections from 5 to 11, the contracting section and spinner 
take the number 12 and 13 respectively, and finally the 
corners from 14 to 17. 

The Figure 2.1 shows the sections and their designated 
numbers, as they are discriminated below: 



SECTION 



DISCRIMINATION 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 



test section 
fan duct 

constant area duct 
constant area duct 



first diffuser (sec 1) 
diffuser (sec 2) 
diffuser (sec 3) 
diffuser (sec 4) 
second diffuser (sec 1) 



contraction cone 
spinner 

first corner with full size vanes 



diffuser (sec 2) 
diffuser (sec 3) 
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Partition of Tunnel Circuit for Evaluation of Losses. 



15 


second 


corner 


with 


full 


size 


vanes 


16 


third 


corner 


with 


full 


size 


vanes 


17 


fourth 


corner 


with 


half 


size 


vanes 



The dimensions of the sections, which are used for 
evaluating the losses, are taken from a WEST COAST RESEARCH 
CO. drawing [Refs. 7,8]. 

In each section the loss K will be a drop in static 
pressure divided by the dynamic pressure at that section, as 
given by RAE & POPE [Ref. 2: p.87] and rewritten below 
Ko = K(q/qo) , 

K = (pi - pf ) / q loss at each section, 

Ko = coefficient of loss related to test section, 
pi = inlet section pressure, 
pf = outlet section pressure, 
q = dynamic pressure at section, 
qo = dynamic pressure at test section. 

In terms of diameter we get 

Ko = K(Do/D) , where 

Do = test section diameter, 

D = local tunnel diameter, 

and finally, for energy ratio of the tunnel, 

ERt = jet energy / summat ion of circuit losses. 

The losses are evaluated by means of a computer program 
written in FORTRAN IV language, which is presented in 
Appendix A, and named LOSS. 

When the evaluation of skin friction is needed, it is 
made by subroutine FRIC in the LOSS program, using the 
equation from RAE & POPE [Ref. 2: p.88], shown immediately 
below, 

1/A = 2 log-^Q Rey A - 0.8, (eqn 2.1) 
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where A. 
number . 


= skin friction 


coefficient and 


Rey 


= Reynolds 


The 


Reynolds number 


at each section 


is 


related to 



Reynolds number at the test section, as follows: 

Rey = VD/v • 

From the continuity equation AV = Ao Vo, and hence 
Rey = (Do/v)Vo (Do/D) , from which is found 

Rey = 2.63E 04 Vo(Do/D) . (eqn 2.2) 

B. CYLINDRICAL SECTIONS 
1 . Test Section 

The losses for this octagonal section are evaluated 
through the eqn 2.3, (eqn 2.44, Ref. 2, p.88) for the 

equivalent cylindrical section, 

Ko= A (L /Be) (Do/De) 4 , (eqn 2.3) 

where De is the diameter of the circle whose circumference 
equals the perimeter of the octagon. 

Then, the dimensions for the test section are: 

Length = 8.00 feet, 

Diameter = 4.36 feet. 

This diameter is the test section diameter (Do) used 
in all the equations related to the test section. 

The losses for the test section at each velocity are 
calculated by the computer program LOSS using the equations 
below, which were obtained from eqn 2. 3, from the equation of 
the Reynolds number, and from the dimensions above. 

section 1: Ko= 1.834A , Rey= 2.63E 04 Vo . 
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2. Fan Duct Losses 



This section begins after the first corner with a 
octagonal cross section, changes to a cylindrical section 
where the two stages operate and ends in another octagonal 
section. From an original wind tunnel drawing (Ref. 8) we 
have the dimensions for this section and, using the same 
procedure as above (test section), we get the following 
dimens ions : 

De = 7.50 ft L= 9.50 ft . 

The same equations used for test section are used 
here and the equations for the computer program are: 

section 2: Ko= . 1448A , Rey = 1.54E 04 Vo . 

3 . Constant Area Ducts 

The losses are evaluated in the same manner as for 
the test section . The dimensions for these two ducts are: 

section 3: De= 13.37 ft, L= 6.00 ft. 

section 4: De= 12.73 ft, L= 6.34 ft. 

The equations for the computer program are: 

section 3: Ko= .0051A , Rey= .86E 04 Vo. 

section 4: Ko= .0067A , Rey= . 90E 04 Vo. 

C. DIFFUSERS 

The diffusers are broken down into parts, the first 
(between test section and first corner) into four, and the 
second (between the second and third corners) into three, 
the purpose of which is to make the evaluation of losses 
more accurately. The diffuser divergence angle, which 
strongly affects the evaluation, is taken for each part as 
the difference between the equivalent small and large 
diameters, divided by the length. 

The equation used comes from RAE & POPE [Ref. 2: p.89], 
and it is 
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Ko= A (X/8tan(a/2) + . 6 tan («/ 2 ) ) , 



( eqn 2.4) 



where A = ( 1- (D1/D2 ) 4 ) (Do/Dl) 4 , and 

A = skin friction coefficient for Reynolds 
number given by eqn 2.2, and based 
on the mean diameter, 



a = divergence angle between walls, 
Do = test section diameter, 

Dl = smaller diameter of diffuser, 

D2 = larger diameter of diffuser. 



1 . First Diffuser 

The diffuser cross section from the test section to 
the first corner varies along its 20 feet of length, and its 
height H and width W are given as a function of the distance 
X (ft), along the flow axis, by the equations: 



H = 3 .52 + .00440 X 2 ' 4 (ft) , 



( eqn 2.5) 



W = 5.02 + .00036 X 3 (ft) . (eqn 2.6) 

Furthermore, its cross-sectional area is: 

Area= .828 W H (ft 2 ) . 

The following table gives the local values of H, W, 
Area and De as function of X for this diffuser, where 

De = equivalent local section diameter (ft). 

The evaluation of the diffuser angle is made by the 
following equations: 

-For the divergence angle of the equivalent conical 
diffuser 



-For the divergence angle between the walls in the 
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TABLE 


I 






PARAMETERS FOR 


DIFFUSER 1 




X 


H 


W 


Area 


De 


0 


3.52 


5.02 


14.63 


4.32 


5.0 


3.73 


5.06 


15 . 64 


4.46 


10.0 


4.62 


5.38 


20.60 


5.12 


15.0 


6.44 


6.28 


33.27 


6.51 


18.7 


8.50 


7.38 


51.94 


8 . 13 



a = 2arc tan( (D2-D1) /2L) . 



( eqn 2.7) 



vertical and horizontal planes respectively: 



=arc tan( (h2-hl) /2L) , 



(eqn 2.8) 



a w =2arc tan ( (w2-wl) /2L) , (eqn 2.9) 

where indexes 1 and 2 refer to small and large local 
sections and L is the length between them. 

For the first diffuser, the equations 2.8 and 2.9 
are used because the angles between the walls are different, 
as shown in Table II . 

For a conservative analysis, the larger values are 
chosen in computing the loss coefficient Ko . 

With the dimensions from this table and the 
equations 2.4 and 2.2 we got the following equations to use 
with the computer program LOSS : 
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TABLE II 

PARAMETERS FOR EACH PART OF DIFFUSER 1 



Part 


L(ft ) 


Dl(ft) 


D2(ft) 


a w (deg. ) 


a h (deg. ) 


1 


5.0 


4.36 


4.46 


0.52 


2.29 


2 


5.0 


4.46 


5 . 12 


3.60 


10.23 


3 


5.0 


5 . 12 


6.51 


9.77 


20.63 


4 


3 . 7 


6.51 


8 . 13 


17.50 


30.85 



section 5 
section 6 
section 7 
section 8 



Ko=.5367A + .00103, 
Ko= . 5416A + .0208 , 
Ko= . 2232A + .0355 , 
Ko= . 0540A + .0196 , 



Rey= 2.60E 04 Vo. 
Rey= 2.40E 04 Vo. 
Rey = 1. 97E 04 Vo. 
Rey= 1.57E 04 Vo. 



2 . Second Diffuser 

For the second diffuser, the values of D1 and D2 are 
substituted into equation 2.7 for evaluation of the 
divergence angle, because the angles between the walls in 
the vertical and horizontal planes have no large difference. 





PARAMETERS 


TABLE 
FOR EACH 


III 

PART OF 


DIFFUSER 2 


Part 


L 


D1 


D2 


a 


1 


9 . 97 


7.40 


7.88 


2.76 


2 


9 . 97 


7.88 


8.89 


5 . 79 


3 


24.10 


8.89 


14.01 


12. 10 



With the dimensions on Table III and equations 2.4 
and 2.2, we got the following equations to use with the 
computer program LOSS. 
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section 9 : Ko= . 14A. + .0004 
section 10: Ko=.09A + .0011 
section 11: Ko=.05A + .0030 

D. CONTRACTION CONE 



Rey= 1.50E 04 Vo. 
Rey= 1.37E 04 Vo. 
Rey = 1.00E 04 Vo. 



This section takes the number 12; it is located between 
corner four and the test section. The equation used is from 
RAE & POPE [Ref. 2: p.91], and is 



Ko= . 32A Lc/Do, 



(eqn 2.10) 



where A = skin friction coefficient for the Reynolds number 
given by eqn 2.2 with D=(Dl + D2)/2; Lc = length of 
contraction cone. 

The equations for the computer program are: 
section 12: Ko=.7338A , Rey= 1.34E 04 Vo. 

E. SPINNER 

This section takes the number 13, and it is located in 
the fan duct. The dimensions are shown in fig 2.2 and the 
loss is evaluated using a method given by Nicolai [Ref. 10: 
p.11-21]. The following equations are for use with the 

computer program: 

Rel = 2.45E-04 Vo, 

Cdf = .455/ (log 10 Rel) , 

Cdo = Cdf + .000616 / Cdf, 

Ko = .00837 Cdo, 

Cdo = drag coefficient , 

Cdf = friction drag coefficient, 

where Rel = Reynolds Number related to body diameter. 

The dimensions are : 

L = Length = 12.0 ft, 

D = Body diameter = 2.0 ft, 

Db = Base body diameter = 0.554 ft. 
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Figure 2.2 Spinner Dimensions. 



F. CORNERS 

The corners are evaluated in three different ways: 

(i) RAE & POPE [Ref. 2 : . pp . 89 , 90] . 

(ii) Bradshaw and Pankhurst [Ref. 5: p.29 ]. 

(iii) Pankhurst and Holder [Ref. 3: p.90]. 

By RAE & POPE the loss equation is 

Ko= (0. 1+ (4.55/ (log 10 Rn) 2 - 58 ) (Do/D) 4 , (eqn 2.11) 



where Rn = Reynolds number based on chord vane, 

D = equivalent diameter of inlet cross section. 

The dimensions are in Table IV, 

With the dimensions from Table IV and the eqn 2.11, the 
equations to be used in the computer program are found and 
are as follow: 



section 14: Ko= 

Rey = 

section 15 : Ko= 

Rey = 

section 16 : Ko= 

Rey = 

section 17 : Ko= 

Rey = 



.00690 + .3122/ (log 10 Rn) 2 - 58 
1.646E 03 Vo . 

.00905+ .4118/ ( log 10 Rn) 2 ' 58 
1.890E 03 Vo . 

.00094+ . 0428/(log 10 Rn) 2 ' 58 
. 609E 03 Vo . 

.00113+ .0515/ (log L0 Rn) 2 ' 58 
. 3 3 4 E 03 Vo. 
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9 
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TABLE IV 

PARAMETERS OF CORNERS 



Section 


Chord ( f t ) 


D ( f t ) 


14 


1.042 


8.52 


13 


1.042 


7 . 95 


16 


1.042 


14.00 


17 


0.521 


13.37 



From Bradshaw and Pankhurst, an empirical equation based 
on tests at Reynolds number 2.0E+05 and 1.9E+06, gives the 
loss coefficient as 



Cp = 1.2 (Uc/v)"' 25 



( eqn 2.12) 



where U = local flow velocity, 
c = vane chord, 
v- kinematic viscosity. 

This equation, when related to the test section velocity 



by means of equation 


Ko= Cp (Do/D)^ yields the 


following 


equations for 


use with 


the computer 


program . 




section 14: 


Ko = 


.0130 Vo"- 25 


> 




section 15 : 


Ko = 


.0160 Vo" • 25 


? 




section 16 : 


Ko = 


.0023 Vo"- 25 


> 




section 17 : 


Ko = 


.0032 Vo" ' 25 


. 




Pankhurst 


and Holder's method 


is a graphic one, which 


gives the variation of loss coefficient with the 


gap / chord 


ratio of the 


vanes . 


The tunnel we are dealing 


with has 


corner - vane 


sections similar to 


type b of figure 40 of 


reference 3 . 


Taking t 


he dimensions 


of the corners 


and vanes 



from the drawing, evaluating the gap/chord ratio for each 
corner and entering the graph on page 93 of reference 3 , 
the value of K ( C is used in Pankhurst & Holder) is found, 
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and used with the equations below: 



K= 2 V H/ p U 2 , ( eqn 2.13) 

Ko=K (Do/D) 4 , (eqn 2.14) 

where U = local velocity at entry of corner, 

VH = drop of pressure across the corner, 

D = equivalent corner diameter at entry, 

K = loss coefficient from the graph. 

Sample calculations 

For corner # 1 we have 
gap / chord = (8.36 - 2.48)/ 12.5'= 0.4704, 

(Do/D) 4 = 0.0967, 

K = 0.235, 

Ko = 0.0227. 

The losses in all four corners, as given by Pankhurst 
and Holder's method are shown in Table V . 

G. SUMMARY OF LOSSES 

The losses evaluated by the computer program for 

sections 1 through 17 for a test section velocity of 200 
Knots, are shown in Table VI ; values for velocities from 

100 to 200 Knots are shown in Table VII . 

The energy ratio given by Ert is multiplied by .9 in 
order to take into account the losses due to leaks and 
joints . 

The results presented in Table VII include, for purposes 
of comparison, the losses evaluated by the other tw ? o 

methods. At the design airspeed of 200 Knots, Pankhurst & 
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TABLE V 

CORNER LOSSES BY PANKHURST AND HOLDER 



corner 


gap / chord 


(Do/D) 4 


K 


Ko 


1 


0.4704 


0.0967 


0.23 


0.023 


2 


0.4256 


0.1334 


0.27 


0.036 


3 


0.4752 


0.0094 


0.23 


0.002 


4 


0.4656 


0.0113 


0.24 


0.003 



Holder's method predicts that 245 HP will be required, while 
Bradshaw & Pankhurst (a later work) indicates a need for 170 
HP. Rae & Pope's method requires 196 HP. The power input to 
the tunnel is 300 x fan ef f iciency ( . 85^ ) = 217 HP. 

The description of the computer program LOSS is 
presented in Appendix D. 
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TABLE VI 

LOSS FOR EACH SECTION OF TUNNEL AT V0=200 KNOTS 



SEC 


ReylO' 6 


A 


Ko 


%Loss 


1 


test 


8.88 


0.0082 


0.0151 


10.99 


2 


cy 1 


5 . 16 


0.0090 


0.0013 


0.95 


3 


cy 1 


2.90 


0.0100 


0.0001 


0.04 


4 


cy 1 


3.04 


0.0100 


0.0001 


0.05 


5 


dif 


8 . 77 


0.0083 


0.0055 


3 . 99 


6 


dif 


8.07 


0.0084 


0.0253 


18.42 


7 


dif 


6 . 65 


0.0090 


0.0375 


27.25 


8 


dif 


5.29 


0.0090 


0.0201 


14.59 


9 


dif 


5.06 


0.0090 


0.0017 


1.21 


10 


dif 


4.62 


0.0091 


0.0019 


1.39 


11 


dif 


3 .37 


0.0100 


0.0035 


2.58 


12 


con 


4.53 


0.0091 


0.0067 


4.86 


13 


spi 






0.0008 


0.60 


14 


cor 


4.56 




0.0067 


4.88 


15 


cor 


4.86 




0.0092 


6 . 70 


16 


cor 


2.77 




0.0009 


0.69 


17 


cor 


2.90 




0.0011 


0.82 
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TABLE VII 

TUNNEL LOSS BY THREE METHODS 



VELOCITY 


TOTAL 


LOSSES 




VELOCITY 

(KNOTS) 


LOSS (B + P) 
(HP) 


LOSS (R + P ) 
(HP) 


LOSS (P+H) 
(HP)' 


100.0 e 


22.2610 


25.4376 


31.2752 


120.0 


38.0835 


43.5826 


53 . 7830 


140.0 


60.1842 


68.9301 


85 . 2725 


160.0 


89.2441 


102 .3171 


126 .8912 


180.0 


125.2569 


143.8930 


179 . 1000 


200.0 


170.4057 


195 .9971 


244 . 5512 
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III. BLADE PERFORMANCE 



A. INTRODUCTION 

The "academic" wind tunnel has two stage counter 
rotating fans, each fed by an electric motor of 150 hp and 
with a fixed speed of 1200 rpm. The main steps used in the 
design procedure were found in the Larson report [Ref. 4,] 
which contains the original blade design. The equations and 
the assumptions made during the design indicate that three 
dimensional effects were not accounted for. For this reason 
the method given by Borst, which is based on the blade 
element approach and cascade theory for determining these 
three dimensional effects, will be used to predict the 
performance of the blades. Usefully, this method requires 
data of two dimensional airfoils only, in order to determine 
the force on each section of the blade. 

So, the blade system with all assumptions made during 
the design, (see section B - Larson report), and the Borst 
method, (see section C ), are the material necessary to 
evaluate the performance of the fans. Knowing the lift 
developed at each blade station (eqn. 3.11) one can find the 
torque (eqn. 3.12) at each station. A blade station is 
defined non-dimensionally by the ratio of the radius of that 
station to the blade tip radius; the root station for this 
tunnel is at .267 and the tip station is at unity. Then 
adding the torques at all stations, and multiplying by the 
rotational speed, one gets the power for the blade. 

B. THE LARSON REPORT 

This report describes the procedure that was followed to 
design the blades, together with the assumptions made during 
the calculations. 

The most important assumptions were: 

- 85% fan efficiency for each stage. 
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of the flow 



Axial velocity distribution 
entering the first stage fan is 



Cm= . 12 Cmr (rt-r)+Cmt , 



(eqn 3.1) 



where Cmr = axial velocity at blade root, 

Cmt = axial velocity at blade tip, 
rt = tip radius , 
r = local radius . 

Maximum test section velocity equal to 200 
knots . 

- Design lift coefficient equal to 0.64, leading 

to the requirement that Cl= .8 . 

- Airfoil type NACA 16 - X12 . 

This design was based on forced vortex flow, in other 
words, the difference in static pressure across the vortex 
is equal to the difference in tangential velocity (dynamic) 
pressure. This leads us to the loading parameter (solidity x 
lift coefficient) and the mean velocity angle, as one can 
see in Appendix A, whose equations are: 



Using these equations and the assumptions mentioned 
above, the dimensions of the blades were calculated, and are 
given on the drawing of the blades (Ref. 9). The blades were 
divided into 11 (eleven) stations from root to tip; for each 
station the chord, thickness, leading edge angle, etc... 



a Cl =( 2 (A-B ) C )/( l+(Cd/Cl) D, and 



(eqn 3.2) 



j3 m = (jSl + jS2 ) / 2 , where 



(eqn 3.3) 
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were determined. 



C. THE BORST METHOD 

This method is based on the blade element approach and 
the vortex theory for determining the three-dimensional 
effects, so that two-dimensional airfoil data can be used 
for determining the resultant force on each blade element. 
By the momentum flow theory, Borst obtained the loading 
parameter (solidity x lift coefficient), as a function of 
the angles j3^, ($2 and 0 m ? that are dependent upon an 
induced angle of attack (a-j_) at values of constant inflow 
angle . 

The loading parameter is : 

ff Cl=( 2 A 2 B (C -D) ) / E, (eqn 3.4) 



where A = 

B = 
C = 
D = 
E = 
Cl = 
Cd = 

01 = 

02 = 
0m = 

r = 

a - 

Now, based on the 



cos 0m ’ 

cosy , 

tanjS! , 

tanj3 2 , 

c°s(0 m -y) , 

lift coefficient, 

drag coefficient, 

apparent inlet angle, 

apparent outlet angle, 

angle of the mean velocity vector, 

atan (Cd/ Cl ) , 

solidity . 

equations 



a i = (jSq ' [$2 )/ 2 ’ an d 



( eqn 3.5) 



0m 01 a i ’ 



(eqn 3.6) 
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and the two dimensional airfoil data and design parameters 
for the blade, one can iterate the equation 3.4 and get the 
induced angle of attack for a specified test section 
velocity and blade station. 

The angle of attack of the blade at any given station 

is : 

a = 01 -j8 + 9 , (eqn 3.7) 

where stagger angle (Ref. 9, Ref. 4 , p.7), which is the 
angle between the chord line of the blade station and the 
line parallel to the axis of rotation of blades, (see fig. 
3.1 ) and 

9- pitch angle (adjustable). 

The Figure 3.1 shows the relation between the variables 
given above. 

D. TORQUE AND POWER EVALUATION 

We will use the Borst method (Ref. 1), to get the 
induced angle for correction of the two dimensional angle of 
attack, and the Larson Report (Ref. 4) to obtain all the 
dimensions of the blades and the preliminary design 
conditions . 

With this material, we are able to follow the procedure 
outlined below to get to the value of torque and power for 
each fan stage. This procedure will lead us to all 
velocities and angles of each blade, finalizing with the 
value of the torque on each blade section. After that, we 
evaluate the torque for the whole blade, and finally the 
torque and power for the fan stage. 

PROCEDURE FOR EVALUATION OF TORQUE AND POWER 

STEP 1) From eqn. 3.1 for a given test section velocity 
and n= 1,200 rpm, calculate W1 =(Cm^ + U^ )l/2> 
where U= 2i r r n. 

As one can see, W1 = f( r,Vo) since Cm = f( r,Vo). In 
Appendix A the axial flow velocity (Cm) was developed as 
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Figure 3.1 Angle of Attack Relations, 
function of Vo . 

STEP 2) Calculate = arc tan ( U / Cm). 

STEP 3) From the drawing of the blades (Ref. 9), the 
values of choi'd are found and then the solidity for the 
blade section is obtained from 

a = Bc/ 27 rr, (eqn 3.8) 

where B = number of blades, c = station chord and r = 
station radius . 

STEP 4) From the "Aerodynamic Characteristics of 24 NACA 
16 - Series Airfoils” - TN 1546 , (Ref. 12, p. 49, 50) the 
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lift and drag coefficient for each section are available at 
MACH number M = Wl / a , where a = sound velocity at 
temperature of 100 degrees Fahrenheit, and for the design 
lift coefficient of 0.64, which was assumed by Larson (Ref. 
4 , p.5,6). The angle of attack used to enter the Cl vsct 
curve is found from eqn.3.7. Find Cl, Cd and evaluate y = 
arc tan(Cd/Cl ) . 

STEP 5) Assuming an initial value of induced angle, 
evaluate the angles given by equations 3.5 and 3.6. 

STEP 6) Evaluate the right side of equation 3.4 using 
the values found in steps 1,2, and 5; evaluate the left side 
of equation 3.4 using the values found in steps 3 and 4. If 
both sides are equal, the value of induced angle assumed is 
correct, otherwise, iterate. 

STEP 7) The corrected angle of attack, found by using 
the induced angle, is seen in figure 3.1 and evaluated as 
follows : 

a cr = a ~ • (eqn 3.9) 

STEP 8) With corrected angle of attack, recalculate the 
step 4 above and get the corrected lift and drag 
coefficient, and evaluate y with these values. 

STEP 9) The mean velocity vector and the lift for each 
blade station are evaluated using the following equations: 

Wm = Cm / cosj3 m , (eqn 3.10) 



Li = 1/2 p Wm 2 Cl c , 



(eqn 3.11) 



where Li = lift at blade station i per unit length, 
Cl = lift coefficient at blade station, 

Wm = mean velocity vector at blade station, 
c = chord at blade station. 
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STEP 10) By definition, the lift force is perpendicular 
to Wm. For evaluation of the torque it is necessary to get 
the component of the resultant force in the plane of the 
fan, as shown in figure 3.2. Then, 

Ri= Li/cos y , 

Fi = Ri cos (/J m -y ) , 

where y = arc tan (Cd/Cl). With Cd and Cl corrected and 
evaluated at step 8), and Li,Ri,Fi are respectively the lift 
force, the resultant force of the lift and drag, and the 
force in the plane of the fan. 




Figure 3.2 Lift Developed at Each Blade Station. 

STEP ll)The force per unit length, Fi, averaged between 
the value at one radial station and next, times the 
difference in radius between these stations, is the force 
used for torque evaluation. This force times the means 
radius of these two stations is the torque. At this point it 
is necessary to go through step 1 to 10 again, in order to 
get Fi for the next station. Then the torque on a blade 
section between two consecutive stations is: 
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Ti =(F i+1 + F i )SR/2, 



( eqn 3.12) 



where S = ( r i+^ - r^) length between stations, 

R = ( r i+ 1 + r i)/ 2 > 

and i refers to the i 1 '* 1 blade station. 

STEP 12) Adding the section torques along the blade, 
multiplying the result by the number of blades and by the 
rotational speed, the power required is found. This power 
provides the test section velocity which was specified at 
the beginning of the procedure. 

E. INTERACTION BETWEEN STAGES 

For evaluating the torque and power for the second 
stage, the procedure given in the previous section can be 
used, but with some adjustments for the different geometric 
characteristics of the blade; for example, the chord, 
stagger angle, etc. . The tangential velocity used is the 
tangential velocity of the. blade at the given station plus 
the tangential velocity of the outflow of the first stage at 
that station. It should be noted that the axial inflow 
velocity is assumed the same as for the first stage. In 
figure 3.4 this interaction can be seen. 

The second stage should be able to accept the tangential 
velocity introduced in the flow by the first stage and add 
the same amount of power as added by the first stage. The 
flow leaving the second stage then has no rotation. 
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Figure 3.3 Sketch for Torque Evaluation at Each Station. 




Figure 3.4 Axial Flow Velocity Triangles. 
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F. CONCLUSION 



Using the computer program BORST and the procedure 
discussed in this chapter, Tables VIII through XX were 
obtained. For all these tables, a pitch angle of 2.4 
degrees was used. 

The Table VIII shows the values of the velocities from 
root to tip, for Vo= 200 Knots, of the first fan stage. 
These velocities are defined as follows: 

U = tangential blade velocity (ft/sec), 

Cm = axial flow velocity (ft/sec), 

W1 = apparent flow velocity (ft/sec). 

The chord at each station is given in feet. 

The Table IX shows the Reynolds and Mach numbers at each 
station, for the first fan stage and at V0 = 200 Knots. It 
is clear that compressibility effects appear at the outer 
20% of the blade only. The maximum Mach number was less than 
.45 and the values of Cl vs Alpha and Cl vs Cd were obtained 
by linear interpolation between the graphs for M= . 3 and 
M= .45 . 

The Table X shows the angles between the axial flow 7 
velocity vector and the velocities relative to the blade, 
the stagger angle, the induced angle, the angle of attack, 
and the lift and drag coefficients. These values are for the 
first fan stage, at each station and at Vo= 200 Knots. The 
values of lift coefficient, are distributed along the 
stations of the blade from .662 to .896, around the uniform 
value of .8, assumed by Larson. These variables are defined 
in Figures 3.1, 3.4, and the angles are given in degrees. As 
expected, at the root, the induced angle is the largest, 
because the greatest interference between the blades occurs 
there . 

The Table XI shows the angle of attack, lift and drag 
coefficients, all with correction for the induced angle. 
Also shown are the lift, the resultant force and the torque. 
This Table refers to each station of the first fan stage, 
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with Vo= 200 Knots. The values of corrected Cl (CLCR) are 
distributed along the blade from .358 to .709, differently 
from that assumed by Larson, which was a uniform .8 . 

The variables in Table XI are defined as follow: 

ALPCR = corrected angle of attack (degrees), 

CLCR = corrected lift coefficient, 

CDCR = corrected drag coefficient, 

DELL = lift force per unit length (lbf/ft), 

DELF = tangential force per unit length (lbf/ft), 

DTOR = torque at middle of two stations (lbf-ft). 

The Tables XII, XIII and XIV, apply to the second fan 
stage. The variables have the same definitions as for the 
first stage. The Table XII shows the inflow velocity 
relative to the blade (W3) larger than Wl. This is due to 
the tangential flow velocity, which has a finite value after 
the first fan stage. At the root the difference between W3 
and Wl is larger than at the tip, this means that the 
tangential flow velocity at the root is larger than at the 
tip. The Table XIII shows the Mach number less than .45, and 
the interpolation described above was again used. The Table 
XIV shows the Cl varying around .8, that is, from .875 to 
.710 as for the first fan stage. 

The Table XV refers to the second stage, with the same 
variables as Table XI plus the variable BETA1-BETA4. This 
Table shows the difference between BETA1 and BETA4 , which is 
around zero. This means that the flow after the second fan 
stage has no tangential velocity. 

The Tables XVI and XVII, are for the uncorrected values 
of the angle of attack, for first stage at Vo=200 Knots. The 
Table XVI shows the values of Cl around .8 . The Table XVII 

shows values of torque (DTOR) larger than those with the 
corrected angle of attack (Table XI). 

The Tables XVIII, XIX and XX are for the second fan 
stage without correction of the angle of attack. A 
comparison of the velocities W3 from Table XVIII, and W3 
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from Table XII shows that the difference, due to the 
correction of the angle of attack has affected the outflow 
of first fan stage. The Table XIX shows the variation of 
BETA3 as a result of the variation of W3 , just mentioned 
above. The variation of Cl was not so large as to affect 
the value of Cd. A comparison of Tables XV and XX shows 
that the interference between the blades reduces the torque 
by about 1/3 at the root and very little at the tip. This is 
further discussed below. 

In order to visualize the effects of the correction on 
the angle of attack, two graphs, each with four curves, were 
plotted; these are: 

1 Figure 3.5; the tangential force on the blade is shewn 
for each station, at a specified test section velocity 
of 200 Knots and at a pitch angle of 2.4 degrees. The 
curves are referred to the first and second stages of 
the blades, with and without correction of the angle of 
attack. This figure shows again the effects of the 
interference. The tangential forces have the largest 
difference at the root, and are practically equal at 
the tip. The data were taken from Tables XI , XV , 
XVII and XX . 

2 Figure 3.6; the angle of attack is shown for each 
station of the blade, at a specified test section 
velocity of 200 Knots and at a pitch angle of 2.4 
degrees. The corrected angle of attack at the root 
became negative and, at the tip, that angle is 
practically equal to the uncorrected one. This again 
shows the strong effects of interference at the blade 
root. The curves are referred as above, and the data 
were taken from Tables X , XI , XIV and XV . 
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Figure 3.5 Tangential Force on the Blade at Each Station 
Vo=200 Knots, Pitch Angle= 2.4 degrees. 



of ATTACK ( DECj 




Figure 3.6 Corrected and Uncorrected Angle of Attack 
Vo=200 Knots, Pitch Angle= 2.4 degrees. 
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TABLE VIII 

VELOCITIES AT EACH STATION, 
FIRST STAGE, V0=200 KNOTS 



STATION CHORD U CM W1 



0 . 


,267 


1. 


. 112 


125 . 


,821 


152. 


606 


197 . 


787 


0 . 


,300 


1. 


,049 


141. 


,372 


150. 


, 340 


206 . 


,369 


0. 


,350 


0. 


, 958 


164. 


,934 


146. 


, 907 


220. 


873 


0. 


,400 


0. 


.875 


188 . 


496 


143 . 


,475 


236 . 


, 887 


0. 


,450 


0, 


.802 


212. 


,058 


140, 


.042 


254. 


, 126 


0. 


, 500 


0. 


.738 


235 . 


. 620 


136 . 


. 609 


272. 


.357 


0. 


.600 


0, 


.631 


282. 


. 744 


129. 


. 743 


311. 


.090 


0. 


.700 


0, 


.546 


329 . 


.868 


122. 


.877 


352. 


.010 


0. 


.800 


0, 


.481 


376 , 


. 992 


116 , 


.011 


394, 


.438 


0, 


. 900 


0, 


.422 


424 , 


. 116 


109 


. 145 


437 , 


. 935 


1 . 


.000 


0 


. 344 


471, 


.240 


102 


.279 


482, 


.212 
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TABLE IX 

REYNOLDS AND MACH NUMBERS AT EACH BLADE STATION, 
FIRST STAGE, V0=200 KNOTS 



STATION 


REY 


MACH 


0.267 


1324694 . 000 


0.171 


0 . 300 


1303483.000 


0. 178 


0.350 


1274545.000 


0.190 


0.400 


1249081.000 


0.204 


0.450 


1227461.000 


0.219 


0.500 


1210513 . 000 


0.235 


0 . 600 


1181581.000 


0.268 


0.700 


1158453.000 


0 . 303 


0.800 


1143157 .000 


0.340 


0.900 


1112776 . 000 


0.378 


1.000 


998700.900 


0.416 
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TABLE XII 

VELOCITIES AT EACH STATION, 
SECOND STAGE, V0=200 KNOTS 



STATION CHORD U CM W3 



0 . 


,267 


0. 


,931 


125 . 


,821 


152. 


, 606 


238 . 


,513 


0 . 


, 300 


0. 


893 


141. 


,372 


150 , 


. 340 


242. 


, 660 


0. 


, 350 


0 . 


,834 


164. 


,934 


146 , 


, 907 


252, 


.430 


0, 


,400 


0, 


, 777 


188 , 


,496 


143, 


,475 


263, 


.640 


0, 


,450 


0. 


. 724 


212, 


,058 


140, 


,042 


276 , 


,489 


0, 


. 500 


0, 


. 673 


235 , 


.620 


136. 


. 609 


292, 


.340 


0 , 


. 600 


0, 


.590 


282, 


. 744 


129 , 


. 743 


327, 


. 947 


0 , 


. 700 


0, 


,519 


329, 


.868 


122, 


,877 


366, 


. 731 


0, 


. 800 


0, 


,461 


376, 


.992 


116, 


.011 


407 


.820 


0 


. 900 


0, 


.411 


424, 


. 116 


109 


. 145 


450 


.087 


1 


.000 


0 , 


.373 


471 


.240 


102 


. 279 


492 


. 649 
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TABLE XIII 



REYNOLDS AND MACH NUMBER AT EACH STATION, 
SECOND STAGE, VO= 200 KNOTS 



STATION 


REY 


MACH 


0.267 


1337683.000 


0.206 


0.300 


1305394.000 


0.209 


0.350 


1268230.000 


0.218 


0.400 


1234027.000 


0.227 


0.450 


1205890.000 


0.238 


0.500 


1185210.000 


0.252 


0.600 


1165001.000 


0.283 


0 . 700 


1146586.000 


0.316 


0.800 


1132560.000 


0.352 


0.900 


1114371.000 


0.388 


1.000 


1106976.000 


0.425 
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TABLE 


XVII 




LIFT, 


RESULTANT FORCE AND 


TORQUE 


WITHOUT CORRECTION OF ANGLE OF ATTACK, FIRST STAGE. 




VO= 200 


KNOTS 




STATION 


DELL 


DELF 


DTOR 


0.267 


38.400 


32.367 


4.133 


0.300 


38 . 362 


30.458 


6.722 


0.350 


39.314 


28.378 


7.101 


0.400 


39.065 


25.482 


7.157 


0.450 


37.845 


22.419 


7.257 


0.500 


39.339 


21.039 


15.625 


0 . 600 


43.683 


19.364 


17 . 160 


0.700 


49.020 


18 . 182 


18.831 


0.800 


56.050 


17.527 


20.424 


0.900 


62.698 


16.647 


20.975 


1.000 


64.297 


14.754 
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TABLE XVIII 

VELOCITIES AT EACH STATION 
WITHOUT CORRECTION OF ANGLE OF ATTACK, SECOND STAGE, 

V0=200 KNOTS 



STATION CHORD U CM W3 



0. 


267 


0. 


931 


125 . 


821 


152. 


606 


232. 


,389 


0. 


.300 


0. 


893 


141. 


,372 


150. 


340 


238 . 


,673 


0 . 


, 350 


0 . 


.834 


164 , 


, 934 


146 . 


907 


250. 


,296 


0. 


.400 


0 . 


,777 


188 . 


,496 


143 . 


,475 


263. 


904 


0. 


,450 


0. 


.724 


212. 


,058 


140. 


, 042 


279. 


,073 


0. 


. 500 


0. 


,673 


235 . 


.620 


136 . 


. 609 


295 . 


.493 


0. 


. 600 


0. 


.590 


282. 


. 744 


129 . 


, 743 


331. 


,213 


0, 


. 700 


0, 


.519 


329. 


.868 


122. 


.877 


369 , 


. 740 


0 , 


. 800 


0. 


.461 


376 , 


. 992 


116 , 


.011 


410, 


.238 


0 , 


. 900 


0 . 


.411 


424. 


. 116 


109, 


. 145 


452, 


. 155 


1 . 


.000 


0, 


.373 


471, 


.240 


102, 


. 279 


495 , 


. 121 
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TABLE XX 

LIFT, RESULTANT FORCE AND TORQUE, SECOND STAGE 
WITHOUT CORRECTION OF ANGLE OF ATTACK, V0=200 KNOTS 



STATION 


DELL 


0.267 


30.579 


0.300 


31.541 


0.350 


31.865 


0.400 


31.182 


0.450 


32.748 


0.500 


36 .370 


0.600 


40.413 


0.700 


46.350 


0.800 


53 .734 


0.900 


61.554 


1.000 


69 . 362 



DELF 
25.789 
25.055 
22.939 
20.289 
19.313 
19.446 
17 . 920 
17 . 190 
16 .797 
16.520 
15 .911 



DTOR 

3.345 
5.484 
5.699 
5.917 
6.472 
14.450 
16.046 
17 . 923 
19.912 
21.663 
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IV. TUNNEL PERFORMANCE 



A. INTRODUCTION 

The tunnel performance depends on losses through the 
tunnel circuit and on the fan efficiency. 

The performance will be evaluated to ascertain how the 
loss around the tunnel circuit is related to the velocity at 
the test section, and if the fans can provide adequate 
energy to generate this test section velocity. The energy 
provided by the blade system depends upon the pitch angle 
which must be limited in order to avoid stall. The pitch 
angle is found for the first and the second stages and 
checked if the flow leaving the second stage has zero 
rotational velocity as required. 

The tunnel performance will be shown through the 
operational envelope of the blade system. 

B. TUNNEL LOSSES 

As we have seen in CHAPTER TWO the losses of each 
section depend on the flow velocity, which implies that they 
depend upon of test section velocity 

Relating the losses for each tunnel section to the test 
section velocity we can get the total losses as a function 
of this velocity, and then relate the latter to flow energy 
required . 

Using the program LOSS with test section velocity 
varying from 100 to 200 Knots, by increments of 20 Knots, we 
get the figure 4.1 which shows the energy required (?jbhp.) 
versus test section velocity (Vo) . 

C. OPERATIONAL ENVELOPE 

Following the method given in Chapter 3, we find the 
power required from the fans as a function of test section 
velocity, with a number of pitch angles of the first and 
second stage blades. 
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Figure 4.1 Energy Required to Overcome Tunnel Losses. 

Since we know the losses around the tunnel in relation 
to test section velocity, we can adjust the pitch angle of 
blades until we get the energy coming into the flow to match 
the losses at same test section velocity. The program BORST 
was used for this. 
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In order to test the sensitivity of the blade flow to 
the approaching velocity profile, we compare the performance 
of the blades in uniform and skewed approach flow. 

Using the axial flow velocity profile, originally used 
by Larson and given by equation 4.1 with the program BORST, 
we got the data shown in Table XXI . This Table presents 
the pitch angle for a specified energy and test section 
velocity, the maximum and minimum pitch angle at that test 
section velocity and its related energy. These values were 
plotted and are shown in Figure 4.2 . This Figure has three 
limit curves, root stall, tip stall and fan efficiency. The 
root stall curve limit is the locus of points at which, at a 
specified test section velocity and pitch angle, the blade 
stalls at the root. The tip stall curve limit is for stall 
at the tip. The fan efficiency curve limit is the locus 
where the operation of the fans is limited by the input 
power and the (assumed) fan efficiency, that is, 300 hp x 

O 

.85 = 217 hp . This means that the power provided by the 
fan system can not exceed 217 hp approximately. The pitch 
angles for both stages are equal at the design speed of 
Vo=200 Knots. 

Using the uniform flow velocity given by equation 4.2 , 
we found the data shown in Table XXII . This Table shows 
that the fans can operate without stall, only at a section 
test velocity from 160 to 180 Knots. These values were 
plotted in Figure 4.3 . This Figure has only the curves of 
root and tip stall. The fans stall before reaching the 
maximum power allowed by the fan efficiency. The blades of 
second fan stage are practically stalled at the root. 

The axial velocity distribution specified for Table 
XXI and Table XXII are developed in Appendix A and shown 
below : 
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Cin = .2034 ( 1 - X ) + .303 Vo , and (eqn 4.1) 



Cm = .363 Vo. 



( eqn 4.2) 



D. FINAL CONCLUSION 

The large subsonic wind tunnel at NPS, powered by two 
counter-rotating fans, never achieved the design 
specif ications . The reason for the poor performance of the 
tunnel was unknown but believed to be due to either poorly 
designed fan blades or to separation in the diffuser. Being 
easier to analyse, the fan blades were chosen for initial 
study. However, to do this analysis, it was necessary to 
compute the losses in the whole tunnel circuit. 

The results can be summarized as follows: 

1) Most of the losses, about 60% of the total, was in the 
first diffuser. The exit end of this diffuser has a high 
value of divergence angle in the vertical plane, which can 
be the cause of a possible flow separation. 

2) When the inflow velocity profile to the fans is assumed 
uniform, the no-stall operational envelope is substantially 
reduced, so there is a much greater likelihood of blade 
stall . 

3) The variation of the pitch angle affects the 
performance of the fans. The pitch angles for both stages 
should be practically the same in order to have no flow 
rotation after the outlet of the second stage, provided the 
inflow velocity profile is the same as that assumed in the 
design. 

4) In conclusion, although this new blade element method 
used to analyse the fan flow predicts some deviation from 
the original isolated blade analysis, it was found that the 
original design was adequate and, therefore, the tunnel 
problem is most probably due to separation in the diffuser. 
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Pitch 

( aec,n f e) 




F igure 4 . 2 



Operational Envelope for Equation 
4 . 1 . 
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Fitch rh fyc - ( te 




Figure 4.3 Operational Envelope for Equation 

4.2. 
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APPENDIX A 



SPECIFIED FUNCTIONAL FORMS 

CORNER 

The evaluation of the corner loss using the turning vane 
loss function presented in appendix B of reference 6 , 

(p.35) with the flow turning angle equal to 90 degrees, 
leads to the equation given by RAE & POPE [Ref. 2: p.90 ]. 

The equation for the turning vane loss coefficient is 

KTV=-1. 605 6 70E- 01+ 1.446 75 3E-02<£ - 2 . 570748E- 04(~> 2 + 2 . 066207 

E-06<£ 3 -6.335764E-09<£ 4 , 

with<£= 90 degrees is found KTV= .15. 

The equation for turning vane loss function is 
K = KTV (2/3+ (1/ 3 ( log 1Q Rnr /log 10 Rn) 2 ' 58 ), 
where 

Rnr = 500,000 (reference Reynolds number, 
according to ref. 2 p.89 ), 
and log-^Q 500,000 = 5.7 . 

Substituting these values in the equation of the turning 
vane loss function, the equation given by Rae & Pope appears 
as follows : 

K = 0.15 (2/3 + 29.7/log 10 Rn) 2 ' 58 ) , 

K = (0.1 + 4.45/ (log 10 Rn) 2 - 58 ), 

Ko = K (Do/D) 4 . 

LOADING PARAMETER (Larson) 

The loading parameter used by Larson in his report comes 
from blade element theory, where the effects of drag are 
considered. This is also shown in Osborne [Ref. 11: p.147], 
by equation 

Vu/Wm= 1/2 Cl( 1+ Cd/Cl cotan (j3-«)) 
where £> = mean blade velocity angle . 

The Larson report (pp. 5,6) gives 

cotan (j3 ~a) = tan (jS m ) , 
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and rearranging the term Vu/Wm as 

Vu/Wm = (U/Cm - (U-Vm)/Cm) Cm/Win , and using the axial 
velocity triangles, the following relations are assured: 
U/Cm= tanj3^, (U-Vu)/Cm = tan ^2 > 

Cm/Wm = cos]3 m 

Substituting these relations into the equation given by 
Osborne, the final equation is 

<jCl = (2 ( tanjS]^- tan ^2 ) cos ]3 m ) / ( 1 + Cd/ Cl tan/3 m ) and is the 
equation used by Larson. 

AXIAL FLOW VELOCITIES. 

Larson assumed in his report the axial flow velocity as: 
Cm = .12 rt (1- r/rt) + Cmt (1) 

where rt= tip radius, r= station radius, Cmt= axial 
velocity at tip blade, and Cmr = axial velocity at the root 
of the blade. Integrating the product of the cross section 
area of the fan duct at the blades times Cm, from root to 
tip, .and setting it equal to the product of the area times 
the velocity of the test section, we find: 

14.9 Vo = 5.465 Cmr + 41.037 Cmt. 

In the equation of Cm, when r is the radius of the blade 
root, Cm = Cmr, and 

Cmr = Cmt / 0.67 . (2) 

Equating (1) and (2), the axial velocity profile for our 
tunnel is found as: 

Cm = ( 0.2034 ( 1- r/rt ) + .303 ) Vo 
UNIFORM AXIAL FLOW VELOCITY. 

This axial profile velocity is used only for purposes of 
comparison. From the continuity equation Cm = Ao Vo / A 
and, taking the known values, we find Cm = .363 Vo 



65 



APPENDIX B 



SAMPLE CALCULATIONS 

Evaluation of the blade station torque. 

For first stage 

Following the 'torque and power' procedure given in 
chapter three, at a specified station the data are as 
follows : 

General data: 

- curves of Cl vsa and Cl vs Cd for air- 
foil NACA 16-X12 at M= . 3 and M=.45. 

a= 1 , 160 ft/sec 

- p = .0023 slug/ft 3 

- v =.000166 ft 2 / sec 

- rt = 3.75 ft 

- Cl des = - 64 

- Vo = 200 Knots = 337.55 ft/sec 

- B = 4 blades 

- Q = 2.4 deg. 

At specified station: 

- X = .3 

- c = 1.049 ft 

- j3 = 40.37 deg . 

STEP 1) Cm = (. 2034(1-. 3 ) + . 303 )x337 . 55 = 150.34 ft/sec 
U = 2x3 . 1416x. 3x3 .75x1200/60 =141.37 ft/sec 
W1 = ( 150 . 34 2 + 141 . 37 2 ) ■'■^ 2 = 206.37 ft/sec 
STEP 2) = arc tan( 141 . 37/ 150 . 34 ) = 43.24 deg. 

STEP 3) o =(4x1.049)/ (2x3 . 14 16x. 3x3.75) = .5936 
STEP 4) M= 206.37/1,160= .18 

a =43.24 - 40.37 + 2.4 = 5.27 deg. 

Entering into curve at M= . 3 are found Cl=.87 and 
Cd = .015 . 



y = arc tan (.015/. 87) = .988 
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5) assume = 5.85 deg. 
j3 2 = 43.24 - 2x5.85 = 31.54 deg. 

j3 m = 43.24 - 5.85 = 37.39 deg. 

6) right side = ( (2xcos 2 37 . 39xcos . 988 (tan43 . 24- 

-tan31.54))/cos(37.39-.988) = 

= .51225 

left side = .5936x.87 =.5164 

The induced angle assumed is adequate and equal 
to = 5.85 deg . 

7) a cr = 5.27 - 5.85 = -.58 deg. 

8) Cl cr =.41 Cd cr =.01 

y cr =arc tan (.01/. 41) = 1.397 deg. 

9) Wm= 150. 34/cos37 . 39 = 189.22 ft/sec 
Li = l/2x. 0023x189. 22 2 x. 41x1. 049=17. 71 lbf/ft 

10) Ri = 17.71 / cosl.397 = 17.71 lbf/ft 
Fi= 17.71 cos (37.39 - 1.397) = 

= 14.33 lbf/ft 
NEXT STATION 
At specified station 
X= .35 
c= .958 ft 

j3 = 46.67 deg. 

STEP 1) Cm = ( .2034(1- .35)+ .303)x337. 55 = 146.90 ft/sec 
U = 2x3 . I4l6x. 35x3 . 75x1200/60 =164.93 ft/sec 
W1 = (146.90 2 +164.93 2 ) 1 / 2 = 220.86 ft/sec 
STEP 2) j3 1 = arc tan( 164 . 93 / 146 . 90 ) = 48.31 deg. 

STEP 3) o =(4x .958)/(2x3.1416x.35x3.75) = .4647 
STEP 4) M= 220.86/1,160= .19 

a=48 . 3 1 - 46.67 + 2.4 = 4.04 deg. 

Entering into curve at M= . 3 are found Cl=.33 and 
Cd = .015 . 

y= arc tan (.015/. 83) =1.033 deg. 



STEP 

STEP 

STEP 

STEP 

STEP 

STEP 
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5) assume a^ = 4.00 deg. 

(^ 2 = 48.31 - 2x4.00 = 40.31 deg. 

£ 3 = 48.31 - 4.00 = 44.31 deg. 

6 ) right side = ( ( 2xcos 2 44 . 3 lxcos 1 . 033 ( t an48 . 3 1- 

-tan40.31))/cos(44. 31-1. 033) = 

= .3864 

left side = .4647x.83 =.3866 

The induced angle assumed is adequate and equal 
to «£= 4.00 deg. 

7) a cr = 4.04 - 4.00 = .04 deg. 

8 ) Cl cr =.455 Cd cr =.01 

y cr =arc tan (.01/. 455) = 1.259 deg. 

9) Wm= 146 . 9 / cos44 . 3 = 205.25 ft/sec 
Li = l/2x. 0023x205. 25 2 x.455x. 958=21. 11 Ibf/ft 

10) Ri = 21.11 / cos 1 . 259 = 21.11 lbf/ft 
Fi= 21.11 cos (44.31 - 1.259)= 

= 15.43 lbf/ft 
For both station at X = .3 and X = .35 
STEP 11) T= (14. 33+15. 43)/2 x3.75x.05x .325x3.75 = 

= 3.40 Ibf-ft 

After that, adding all station torques and multiplying 
by the rotational speed, the power is found. 

For the second stage, the same procedure can be 
followed with some adjustments as mentioned before in 
chapter 3, section E. 



STEP 

STEP 

STEP 

STEP 

STEP 

STEP 
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APPENDIX C 



DESCRIPTION OF COMPUTER PROGRAMS 

The following pages contain the FORTRAN PROGRAMS 
developed to evaluate the blade performance and the tunnel 
losses . 

The LOSS FORTRAN is a program to evaluate the losses 
around the circuit of the large Academic Wind Tunnel of NPS . 
This program uses one subroutine called SKIN to calculate 
the skin friction when the Reynolds number is given; it does 
this by means of equation 2.1. The main program uses the 
equations for Ko and Rey for each section of the tunnel, 
given in chapter two. The meaning of variables used in the 
program are given at beginning of the main program and 
subroutine. The listing of the program is presented in 
Appendix D. 

The BORST FORTRAN is a program used to evaluate the 
performance of the blades of both stages. This program uses 
the procedure presented in chapter three section D, to 
evaluate the torque and power of the first stage, and 
section E to evaluate the torque and power of the second ; 
it uses six subroutines as follows: FIRST is a subroutine 
to evaluate the parameters of the blade of each station at a 
specified test section velocity, for the first stage. 
SECOND is a subroutine like FIRST but it evaluates for the 
second stage. CLNEW is a subroutine to compute the lift 
coefficient for each section of the blade using Newton's 
Forward and Backward Interpolation Formula, with the data 
given by curves of Cl vs a in reference 12 . CDNEW is a 
subroutine to compute the drag coefficient for each section 
using the same formula as in CLNEW, with data given by the 
curve of Cd vs Cl in the same reference as in CLNEW. INDUC 
is a subroutine used to compute the induced angle of attack 



69 



through iterations using the equation 3.4. LIFT is a 
subroutine to compute the torque and power for second stage 
of blades calling all the subroutines except FIRST; this 
subroutine is called by the main program after it evaluates 
the torque and power of the first stage. 

This program is listed in Appendix E and a description 
of parameters appears at the beginning of either the main 
program or the subroutine. 
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COMPUTER PROGRAM ’ LOSS 
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RELATED TO CORNER EVALUATION BY PANKHURST & HOLDER 
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BLOCK TO INTERPOLATE POINT IN CD VS CL CURVE USING 
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